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Abstract 

In the approach, based on analyticity and unitarity and assuming an influence of 
coupled channels, experimental data on the isovector P-wave of tttt scattering have been 
^ . analyzed to study /5-like mesons below 1.9 GeV. 
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The investigation of vector mesons is an actual problem up to now due to their role in forming 
. the electromagnetic structure of particles and because, e.g., in the p-family, only the p(770) 
meson can be deemed to be well understood [1]. The other p-like mesons must be either still 
\ confirmed in various experiments and analyses or their parameters essentially corrected. For 
example, the p(1250) meson was discussed actively some time ago [2, 3] and it was confirmed 
relatively recently in the amplitude analysis of the LASS Collaboration [4] and in combined 
analysis of several processes [5]. However this state is referred to only slightly in the PDG 
issue [1] (the relevant observations are listed under the p(1450)). 

On the other hand, the irn interaction plays a central role in physics of strongly interacting 
particles and, therefore, it has always been an object of continuous investigation. Let us note 
only some recent works devoted also to the theoretical study of the isovector P-wave of txtx 
■ scattering. First, there are the analyses of available experimental data on the nn scattering 
>• , utilizing the Roy equations [6, 7, 8] and the forward dispersion relations [9, 10], in which, e.g., 
the low-energy parameters of the tttt scattering were obtained. Second, there are the works in 
which the low-energy parameters are calculated in chiral theories with the linear realization 
of chiral symmetry [11, 12]. 

We have used our model-independent method [13] based on the first principles (analyticity 
and unitarity) directly applied to analysis of experimental data, aiming at studying the p-like 
mesons below 1.9 GeV and obtaining the 7T7r-scattering length. Unfortunately this method, 
using essentially a uniformizing variable, is applicable only in the 2-channel case. Here the 
7T7T and U7T channels are allowed for explicitly (in the threshold range of the latter, one has 
observed a deviation from elasticity of the P-wave tttt scattering). Influence of other coupled 
channels is supposed to be taken into account through the background. In order to investigate 
the coupling of resonances with these other channels, we apply also multichannel Breit-Wigner 
forms to generate the resonance poles. 

The paper is organized as follows. In Section II, we outline the method of the uniformizing 
variable in applying it to studying the 2-channel tttt scattering and present results of the 
analysis of the available data [14]-[16] on the isovector P-wave of nn scattering. Section III 
is devoted to analysis of the same data using the Breit-Wigner forms. Finally, in Section VI, 
we summarize and discuss obtained results. 
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2 Method of the uniformizing variable 



Let the 7T7r-scattering S'-matrix be determined on the 4-sheeted Riemann surface with the 
right-hand branch-points at 4m 2 + and (m w + m^o) 2 and also with the left-hand one at s = 0. 
It is supposed that influence of other branch points can be taken into account through the 
background. The Riemann-surface sheets are numbered according to the signs of analytic 
continuations of the channel momenta 



as follows: signs(Im/c 1 , lmk 2 ) = ++, — h, , H — correspond to sheets I, II, III, IV, respec- 
tively. 

The ^-matrix is supposed to be S — S res Sb g where S res represents resonances and Sb g , the 
background. In general, an explicit allowance for the (m w + m 7T o) 2 branch point would permit 
us to describe transitions between the tttt and ujtt initial and final states with the help of 
the only one function d(ki, k 2 ) (the Jost matrix determinant) using the Le Couteur-Newton 
relations [17]. Unfortunately, data on the process tttt — > ujtt are absent. 

In Ref. [13] it was shown how one can obtain the multichannel resonance representations 
by poles and zeros on the Riemann surface with the help of the formulae, expressing analytic 
continuations of the matrix elements, describing the coupled processes, to unphysical sheets 
in terms of those on sheet I. It is convenient to start from resonance zeros on sheet I. Then 
in the 2-channel nn scattering, we have three types of resonances: (a) described by a pair of 
complex conjugate zeros in the S- matrix element on sheet I and by a pair of conjugate shifted 
zeros on sheet IV; (b) described by a pair of conjugate zeros on sheet III and by a pair of 
conjugate shifted zeros on sheet IV; (c) which correspond to a pair of conjugate zeros on sheet 
I, a pair on sheet III and two pairs of conjugate zeros on sheet IV. The poles on sheet II, III 
and IV are situated in the same energy points as the corresponding zeros on sheet I, IV and 
III, respectively. Note that the size of shift of zeros on sheet IV relative to the ones on sheets 
I and III is determined by the strength of coupling of the channels (here tttt and ujtt). The 
cluster kind is related to the nature of resonance. 

With the help of the uniformizing variable 1 



the considered 4-sheeted Riemann surface is mapped onto the f-plane, divided into two parts 
by the unit circle centered at the origin. Sheets I (II), III (IV) are mapped onto the exterior 
(interior) of the unit disk on the upper and lower f-half-plane, respectively. The physical 
region extends from the point % on the imaginary axis (tttc threshold) along the unit circle 
clockwise in the 1st quadrant to the point 1 on the real axis (utt threshold) and then along 



the real axis to the point b = J [m^ + m w o + 2m 7r +)/(m a; + m n o — 2m n +) into which s = oo 



1 The analogous uniformizing variable has been used, e.g., in Ref. [18] in studying the forward elastic pp 
scattering amplitude and in Ref. [19] in the combined analysis of data on processes tttt — > irn, KK in the 
channel with I G J PC = 0+0++. 
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is mapped on the v-plane. The intervals (—00, —b], [— b^ 1 , b^ 1 ], [b, 00) on the real axis are the 
images of the corresponding edges of the left-hand cut of the 7T7r-scattering amplitude. The 
(a) resonance is represented in S(irn — > nn) by two pairs of poles on the images of sheets II 
and III, symmetric to each other with respect to the imaginary axis, and by zeros, symmetric 
to these poles with respect to the unit circle. Note that the symmetry of the zeros and poles 
with respect to the imaginary axis appears due to the real analyticity of the S'-matrix, and 
the symmetry of the poles and zeros with respect to the unit circle ensures a realization of the 
known experimental fact that the ixix interaction is practically elastic up to a vicinity of the 
uj7i° threshold. 

The resonance part of S'-matrix S res becomes a one-valued function on the f-plane and, in 
the 7T7T channel, it is expressed through the (i(v)-function as follows 2 : 

S -t^ (2) 
* res ~ d(v) [ ) 

where d(v) represents the contribution of resonances, described by one of three types of the 
pole clusters in the 2-channel case, i.e., 

M 

d = v~ M H(l - v*v)(l + v n v) (3) 

n=l 

with M the number of pairs of the conjugate zeros. 
The background part Sb g is taken in the form: 



-> bg = exp 



2l I W S 4m ' + I ( ai + a 2 6(s - si) + a 3 0(8 - 8 2 ) 



(4) 



where ctj = + ibi, Si is the threshold of 4n channel noticeable in the p-like meson decays, 
s 2 is the threshold of p2ir channel. Due to allowing for the left-hand branch-point at s = in 
the uniformizing variable (1), a± — b± — 0. Furthermore, b 2 = is an experimental fact. 

With formulas (2)-(4), we have analyzed data [14]— [16] for the inelasticity parameter (77) 
and phase shift of the 7T7r-scattering amplitude (5): S(tttt — > im) = i]exp(2i5). First, we 
considered the data of Refs. [14, 15] introducing three (p(770), p(1250 — 1580) and p(1550 — 
1780)) and four (the indicated ones plus p(1860 — 1910)) resonances. 

From a variety of possible resonance representations by pole-clusters, the analysis selects 
the following one to be the most relevant: the p(770) is described by the cluster of type 
(a) and the other resonances by type (b). Description in both assumed cases is satisfactory: 
X 2 /NDF is 204.795/(146-15) = 1.563 (three resonances) and 197.936/(146-19) = 1.558 (four 
resonances). The background parameters are: a 2 = 0.0039, a 3 = 0.0253, and 63 = —0.0337 
for the case of three resonances and a 2 = —0.00065, a 3 = 0.00432, and 63 = 0.0001 for four 
resonances. When calculating x 2 f° r the inelasticity parameter, three points of data [15] at 
990, 1506 and 1825 MeV have been omitted in both cases as giving the anomalously big 



2 Other authors also have used the parameterizations with the Jost functions at analyzing the S'-wave tttt 
scattering in the one-channel approach [20] and in the two-channel one [21]. In latter work, the uniformizing 
variable hi has been used and the 7T7r-threshold branch-point has been neglected, therefore, their approach 
cannot be employed near by the tttt threshold. 
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contribution to x 2 . When calculating x 2 for the phase shift, three points of data [15] at 1841, 
1804 and 1882 MeV have been omitted in the case of three resonances and those at 765, 1643 
and 1804 MeV in the case of four resonances. 

Let us show the obtained pole clusters on the lower y^-half-plane in Table 1 (it is clear 
that there are also complex-conjugate poles on the upper half-plane). 



Table 1: Pole clusters for the p-like resonances. 



Three resonances 


Sheet 


II 


III 


IV 


p(770) 


y/s;, MeV 


767-i(74.8) 


786.5 - i(64.1) 




p(1250) 


y/s~ r , MeV 




1249.9 -^(156.5) 


1249.6 -^(153.1) 


p(1600) 


MeV 




1576-^(141.7) 


1575.6 -^(78.4) 


Four resonances 


Sheet 


II 


III 


IV 


p(770) 


y/s;, MeV 


767.2 - i(74.8) 


785.2 -^(64.1) 




p(1250) 


y/s;, MeV 




1249.5 -^(156.4) 


1249.7-^(152.2) 


p(1600) 


y/s~ r , MeV 




1578.7-^(142.5) 


1578.2 -^(77.1) 


p(1900) 


y/s;, MeV 




1871.5-i(95.7) 


1894.9-i(98.2) 



On figures 1, we demonstrate results from our fitting to data [14, 15]. The dashed curve is 
for the three-resonance description and the solid one for the four-resonance case. 



77 + 71 -> 77 + 71 



77 + 77 -> 71 + 7T 
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Figure 1: The phase shift of amplitude (the left part) and module of matrix element (the right 
part) of the P-wave 7T7T-scattering. The dashed curve - 3-resonance description; the solid curve 
- 4-resonance one. The data are from Refs. [14, 15]. 

Though x 2 /NDF is practically the same in both cases, consideration of the obtained pa- 
rameters and energy dependence of the fitted quantities suggests that the fourth resonance 
p(1900) is desired. Therefore, in the 4-resonance picture, we have also carried out an analysis 
of data [16] jointly with the already considered data [14, 15]. In Ref.[16], results of two anal- 
yses are cited: one uses the s-channel helicity amplitude when extracting the 7r7T-scattering 
amplitude on the vT-exchange pole; in the other, the t-channel one is used instead. Therefore, 
we have taken both analyses as independent. There are given the data for the phase shift of 
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amplitude below the KK threshold. Comparing these data with the ones of Refs.[14, 15], one 
can see that the points of the former lie systematically by l°-5° higher than the ones of the 
latter, except for two points of Ref.[14] at 710 and 730 MeV, which lie by about 2° higher 
than the corresponding points of Ref. [16] and which are omitted in the subsequent analyses. 
Since we do not know the energy dependence of the remarked deviations of points, we have 
supposed a constant systematic error that must be determined in the combined analysis of 
data. We have obtained a satisfactory description with xV NDF = 284.777/(186-19) = 1.705 
and with the indicated systematic error equal to —1.9°. 

Let us indicate (Table 2) the obtained resonance zero positions on the right-hand f-half- 
plane (there are also zeros symmetric to the indicated ones with respect to the imaginary 
axis) . 



Table 2: Positions of the resonance zeros on the right-hand v-half-plane. 



Resonance 




p(770) 


Vl = 1.044245 + 0.215405i v 2 = 0.926647 - 0.178156i 


p(1250) 


v 3 = 1.238348 - 0.032333i v 4 = 0.806223 - 0.021908i 


p(1600) 


v b = 1.282453 - 0.007478? v 6 = 1.282453 - 0.007478i 


p(1900) 


v 7 = 1.305758 - 0.005171? v$ = 0.766563 - 0.003215i 



The background parameters are: a 2 = 0.00599, 03 = 0.02585 and 63 = —0.00498. 
On figures 2, we demonstrate energy dependences of the analyzed quantities compared 
with all used experimental data [14]-[16]. 
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Figure 2: The phase shift of amplitude and module of matrix element of the P-wave tttt- 
scattering. The data from Refs [14]- [16]. 

In Table 3, we show the obtained pole clusters on the lower ^-half-plane. 

Masses and widths of the obtained p-states can be calculated from the pole positions on 
sheet II for resonances of type (a) and on sheet IV for resonances of type (b) . If the resonance 
part of the amplitude reads as 



Table 3: Pole clusters for the p-like resonances. 



Sheet 


II 


III 


IV 


p(770) 


y/s;, MeV 


767.1 -^(73.3) 


783 - i(65.8) 




p(1250) 


y/T r , MeV 




1251.2-i(152.2) 


1249.5 -i(144.6) 


p(1600) 


y/s^ : MeV 




1584.6-i(141.3) 


1579.7-^(75) 


p(1900) 


y/s;, MeV 




1872.3-i(96.6) 


1895.5 -^(95) 



we obtain for the masses and total widths, respectively, the following values (in the MeV 
units): 



for 


p(770), 


770.6 


and 


146.6; 


for 


p(1250), 


1257.8 


and 


289.2; 


for 


p(1600), 


1581.5 


and 


150; 


for 


p(1900), 


1897.9 


and 


190. 



3 The Breit-Wigner analysis of P-wave nn scattering 

In various works [1], it was shown that the p-like resonances obtained in the previous section 
have also other considerable decay channels in addition to those considered explicitly above. 
It was observed that the p(1450) (and/or p(1250)?) have also the decay modes: rjp° (< 4%), 
47r (seen), and (fin (< 1%). The p(1700) has also the large branching into the 4n (large), p2n 
(dominant), and rjp° (seen) channels. 

To include explicitly influence of several interested channels and to obtain information 
about couplings with these channels only on the basis of analysis of the 7T7r-scattering data, 
we have used 5-channel Breit-Wigner forms in constructing the Jost matrix determinant 
d(ki, • • • , k 5 ), to generate the resonance poles and zeros in the Le Couteur-Newton relation: 

_ d(—ki, ■ ■ ■ ,k 5 ) 
bres ~ d(k ir ..,k 5 ) ■ [b) 

In eq.(5) ki, &2, k^, k± and k 5 are the nn-, 7i + n^2ii -, 2(ti + it~)-, rfln- and cJ7r°-channel mo- 
menta, respectively. The ci-function is taken as d — d res db g with d res , describing resonance 
contributions, and dt> g , the background. 

The Breit-Wigner form for the resonance part of the rf-function is assumed as 



d r es(s) = n 



5 

2 

rj ■* "''j Jr j 

3=1 



M 2 r -S-iY,P% Rrj /, 



(6) 



where p r j = kj(s)/kj(M 2 ), f^/M r is the partial width of resonance with mass M r , Rrj is a 
Blatt-Weisskopf barrier factor [22] conditioned by the resonance spins. For the vector particle 
this factor have the form: 

1 + \Us-Amj r rj y 
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with radius r r j = 0.705 fermi for all resonances in all channels as a result of our analysis. 
Furthermore, we have assumed that the widths of resonance decays to it + it~2-k and 2(7r + 7r~) 
channels are related with each other by relation: f r2 = f r3 /V2. 
The background part db g is 



d, 



bg = exp 



Am 2 



Ot\ + Oi 2 



«1 



0(s - Sl ) 



(8) 



where on = at + ibi, si is the threshold of pin channel (it is clear that a 2 and b 2 take into 
account also influence of other channels opened at higher energies than the pin threshold); b\ 
is taken to be zero. 

With formulas (5)-(8) we have carried out the analysis, just as in the previous section, 
both with three resonances and with the four ones. We have obtained the same reasonable 
description in both cases: the total x 2 /NDF = 316.206/(186 — 17) = 1.871 for three resonances 
and equals 315.254/(186 — 22) = 1.922 for four resonances. Let us show the obtained values 
of resonance parameters for the second case in Table 4. The systematic error of data [16], 
discussed in the previous section, is equal to —2.031° in this case. When calculating x 2 for the 
inelasticity parameter, three points of data [15] at 990, 1030 and 1825 MeV have been omitted 
as giving the anomalously big contribution to \ 2 ■ When calculating x 2 for the phase shift, three 
points of data [16] have been omitted: one at 790 MeV from the s-channel analysis, and two at 
790 and 850 MeV from the t-channel one. For the background we find: a\ = —0.001 ± 0.002, 
a 2 = -0.1065 ± 0.012, and b 2 = 0.00715 ± 0.016. 



Table 4: The p-like resonance parameters (all in the MeV units). 



Resonance 


M 


frl 


frl 


fr3 


fri 


fr5 




p(770) 


777.57±0.33 


343.5±0.75 


25.8±5.9 


36.5±8.4 






«154.3 


p(1250) 


1247.8±15.9 


75.9±7.5 


175.8±40.7 


248.7±57.5 


176±112 


162±116 


>124 


p(1600) 


1580±4.7 


247±5.3 


236.8±8.7 


334.9±12.3 


159±32 


106.5±40.6 


>168 


p(1900) 


1908.3±38.3 


46.4±16 


138 


195 


29.7 


34 


>32 



Note that, in Table 4, we do not show the errors for / 42 , / 43 , / 44 , and / 45 , because insuffi- 
ciency of data does not permit us to fix reliably these parameters. On figure 3, we demonstrate 
results from our fitting to data [14]- [16]. We may conclude that whereas the model-independent 
analysis testifies somehow in favour of existence of the p(1900), the Breit-Wigner approach 
cannot verify this result. 

We have also calculated the isovector P-wave length of tttt scattering, a\. Its value is 
shown in Table 5 in comparison with the ones known from various evaluations in the local [12] 
and non-local [11] Nambu-Jona-Lasinio (NJL) model and from the ones with the use of Roy's 
equations [9, 8, 7] (in the second work one apply also the chiral perturbation theory (ChPT)to 
construct a precise 7T7r-scattering amplitude at s 1 / 2 < 0.8 GeV). 



4 Conclusions 

The reasonable description of all the accessible experimental data on the isovector P-wave 
of 7T7T scattering for the inelasticity parameter (rj) and phase shift of amplitude (5) [14]- [16] 
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7T + 71 -> 7T + 71 



71 + 7T -> 7T H- 7T 




Figure 3: The phase shift of amplitude and module of matrix element of the P-wave tttt- 
scattering. The curves show result of fitting to data [14]- [16] using the Breit-Wigner form. 

Table 5: Comparison of values of the nn scattering length a\ from various approaches. 



a\[10-^m-i] 


References 


Remarks 


33.98 ±2.03 


This paper 


Breit-Wigner analysis 


34 


[12] 


Local NJL model 


37 


[11] 


Non-local NJL model 


37.9 ±0.5 


[8] 


Roy equations using ChPT 


38.4 ±0.8 


[9] 


Roy equations 


39.6 ±2.4 


[7] 


Roy equations 



have been obtained up to 1.88 GeV based on the first principles (analyticity and unitar- 
ity) directly applied to analysis of the data. Analysis has been carried out in the model- 
independent approach using the uniformizing variable (here the satisfactory description is 
obtained: x 2 /NDF = 1.558) and applying multichannel Breit-Wigner forms to generate the 
resonance poles and zeros in the S- matrix (x 2 /NDF = 1.922). The aim of analysis (except for 
obtaining a unified formula for the P-wave nn scattering amplitude in the whole of investi- 
gated energy range) was to study the p-like mesons below 1.9 GeV and to obtain the P-wave 
7T7r-scattering length. 

For the p(770), the obtained value of mass is a little bit smaller in the model-independent 
approach (770.6 MeV) and a little bigger in the Breit-Wigner (777.57 ± 0.33 MeV) one than 
the averaged one (775.5 ± 0.4 MeV) cited in the PDG tables [1], however, it also occurs in 
analysis of some reactions [1]. The obtained value of the total width in the first case (146.6) 
coincides with the averaged PDG one (146.4 ± 1.1 MeV) and it is a little bit bigger in the 
second case (~ 154.3 MeV) than the averaged PDG value, however, this is encountered also 
in other analyses [1]. Note that predicted widths of the p(770) decays to the 47r-modes are 
significantly larger than, e.g., the ones evaluated in the chiral model of some mesons based on 
the hidden local symmetry added with the anomalous terms [23]. 

The second p-like meson has the mass 1257.8 MeV in the 1st analysis and 1247.8 ± 15.9 in 
the 2nd one. This differs significantly from the mass (1459 ±11) of the 2nd p-like meson cited 
in the PDG tables [1]. We told already in Introduction that the p(1250) meson was discussed 
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actively some time ago [2, 3], and next the evidence for it was obtained in some analyses [4, 5]. 
To the point, if this state is interpreted as the first radial excitation of the 1 + 1 state, then 
it lies down well on the corresponding linear trajectory with a universal slope on the (n, M 2 ) 
[24] (n is the radial quantum number of the qq state), while the meson with mass M = 1450 
MeV turns out to be considerably higher than this trajectory. 

The third p-like meson turns out to have the mass 1580 MeV rather than 1720 MeV cited 
in the PDG tables [1]. Note that in a number of previous analyses of some reactions one has 
also obtained the resonance with mass near 1580 MeV [1]. However, some time ago it was 
shown that the 1600- MeV region contains in fact two p-like mesons. This was made on the 
basis of investigation of the consistency of the In and 4-7T electromagnetic form factors and 
the 7T7r-scattering length [25] and as a result of combined analysis of data on the 2n and 4-7T 
final states in the e + e~ annihilation and photoproduction [26]. We assume this possibility, 
i.e., that in the energy range 1200-1800 MeV there are three p-like mesons, but for the final 
conclusion the combined analysis of several processes with which the investigated resonances 
are appreciably coupled have to be performed. Note also a rather big obtained coupling of 
these p-like mesons with the 47r channels. 

Finally, as to the p(1900), in this energy region there are practically no data on the P- 
wave of 7T7T scattering. The model-independent analysis, maybe, somehow testifies in favour 
of existence of this state, whereas the Breit-Wigner one gives the same description with and 
without the p(1900). For more definite conclusion about this state, the P-wave nn scattering 
data above 1.88 GeV are needed. Furthermore, the combined analysis of coupled processes 
should be carried out. 

Yu.S. acknowledges support provided by the Votruba-Blokhintsev Program for Theoretical 
Physics of the Committee for Cooperation of the Czech Republic with JINR, Dubna. P.B. 
thanks the Grant Agency of the Czech Republic, Grant No.202/05/2142. 
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